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Using exact diagonalization we examine the response of several most prominent fractional quan- 
tum Hall states to a single local impurity. The 2/3 singlet state is found to be more inert than 
the polarized one in spite of its smaller incompressibility gap. Based on its spin-spin correlation 
functions we interpret it as a liquid of electron pairs with opposite spin. A comparison of different 
types of impurities, non-magnetic and magnetic, is presented. 
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The effect of impurities on the incompressible ground 
states was in the focus of many theoretical studies since 
the early days of the fractional quantum Hall effect 
(FQHE). The principal interest is based on the knowl- 
edge that the FQHE occurs only in high mobility sam- 
ples implying that random electric potential ubiquitous 
in experiments leads to a decrease of the incompressibil- 
ity gap eventually destroying the incompressibility com- 
pletely. A large part of the theoretical work therefore 
deals with random ensembles of impurities - and this, due 
to the complexity of such calculations, in context of the 
most robust incompressible FQH state, bearing the name 
of R. Laughlin, which occurs at filling factor— v = i. 

A single impurity can, however, be also viewed as a 
probe into the correlated state revealing additional in- 
formation on its properties and nature 3 -. Beyond this, 
as we believe fundamental interest, this information 
may help to understand new correlated states for which 
only numerical many-body wavefunctions are available. 
The particular example we bear in mind are the half- 
polarized states at v = | and | discovered first in optical 
experiments^. Numerically, the wavefunctions of suitable 
candidates can be obtained by exact diagonalization but 
it was not possible to link them conclusively with any 
of the physical states proposed subsequentl y 5 ' 6 ' 78 ^. The 
current study of the parent polarized and unpolarized 
states may provide a guideline to achieve this. 

First, the Laughlin wavefunction is investigated and 
different regimes regarding tis response to a point impu- 
rity are identified. Here, we build on the seminal articles 
of Zhangi^ and Rezayiii. A comparison between point 
and <5-linc impurities, is presented. Proceeding to filling 
factors | and |, it proves advantageous to apply dif- 
ferent types of impurities, magnetic and non-magnetic. 
This investigation reveals basic structural properties of 
the various FQH states like local density and polariza- 
tion responses to impurities. Based on these, we propose 
that the v = | singlet state is a v = 1 liquid of electron 
pairs with opposite spin. 

We recall that the three filling factors studied here, |, 



| and | are all closely related. In the composite fermion 
(CF) picture^, they correspond to filling factors v* equal 
to 1, 2 and 2 respectively and the last two differ only 
in the direction of the effective field. Systems at v = | 
and | can therefore be expected to behave very similarly, 
at least within the mean field of the CF. On the other 
hand, with the spin degree of freedom frozen out (e.g. by 
large Zeeman energy), the ground states at v = ^ and | 
are identical up to a particle-hole conjugatio n 13 ' 14 . Here 
we show how far these statements of close relationship 
apply when the excitations become involved in addition 
to the ground states as it is the case with inhomogeneous 
systems. 



I. THE MODEL AND POINT IMPURITIES 

It is a common notion that the Laughlin state is in- 
compressible. This statement however relates to a ther- 
modynamical property and does not contradict the fact 
that even an arbitrarily small impurity potential locally 
changes the electron density, Fig. Q] 

As a function of the impurity strength Vq, the local 
density at the position of the impurity n(ro) shows three 
distinct types of behaviour marked by I, II and III in 
Fig. [21 The boundary between I and II can roughly be 
identified with the incompressibility gap E g . While the 
density change An(ro) is roughly proportional to Vq for 
| Vol ^ Eg (region II), reminiscent of a standard compress- 
ible behaviour as of a Fermi gas or liquid, the density 
change is proportional to Vq for weak impurities (region 
I) . The latter non-linear region does not appear in earlier 
datai 1 - on a sphere and we attribute it to the center-of- 
mass part of the wavefunction as we explain below. Fi- 
nally, the response diverts again from the linear regime 
for very strong impurities (region III). We will not inves- 
tigate this regime here and focus only on the regions of 
An(r ) oc V and |An(r )| <x V 2 . 

As a tool of study we use the exact diagonaliza- 
tion (ED) with Coulomb-interacting electrons on a 



2 




FIG. 1: Right: electron density n(x, y) of the Laughlin state 
in a square with periodic boundary conditions under the in- 
fluence of a single point impurity of intermediate strength. 
Six electrons (implying a = b = 10. 6^0 to ensure v = |), 
Vo = 0.3e 2 /(47re4), a = 1.0f (Eq. H}. Left: section of 
n(x,y) along the diagonal (ED), comparison to the simple 
model based on x{<t) = 8(\q — qo\), Qo^o = 1.4, see text. 
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FIG. 2: Compression of the Laughlin state by a point impurity 
measured by the local density at vq (Eq. [lj. Strength of the 
impurity Vo is divided into regions I, II and III described in 
text. 

toru o 2 ' 15 i 16 . Contrary to the spherical geometry, a point 
impurity 

V(r) = V Q exp [-(?- r ) 2 /a 2 ] (1) 

does break geometric symmetries on the torus described 
by the quantum numbers^ k = (k x , k y ). This makes the 
study more difficult from the computational view (large 
dimension of the Hilbert space) but it also gives the sys- 
tem the full freedom of choosing the ground state. The 
density response n(x, y), Fig. [T], clearly follows the impu- 
rity form (rotationally symmetric) on short ranges and 
it is deformed by the periodic boundary conditions on 
distances comparable to the size of the elementary cell 2 - 
a = b = \f2-KN m £(). Here £q is the magnetic length and 
N e , N m are the number of electrons and the number of 
magnetic flux quanta. 

It should be noted that the oscillations in the density 



response n(x, y) arc not of the Fricdcl type known in a 
Fermi ga o 18 i 19 ' 20 which occur as interferences at the edge 
of the Fermi-Dirac distribution at the Fermi wave vector. 
In the regime of linear response, the oscillations like in 
Fig. [2] can be described up to a very good precision by 
a model dielectric response function x(o) = ~ Qo)- 
This is an approximation based on the observation that 
the magnetoroton minimum at qo£o ~ 1-4 is the lowest 
no-spinffip excitation at v = ^. A more realistic x(o) 
based on the single mode approximation was given by 
MacDonald3i and we recall that it is very different from 
x(q) of a Fermi ga a 18 ' 21 . 

It was noted already by Rezayi and Haldaneii that 
the form of the density response calculated for the lin- 
ear regime, An(r) cx Jo(rqo), remains almost unchanged 
even in the non- linear regime |Vo| 2> E g . This state- 
ment applies also for point impurities on a torus with 
two comments, (a) The density profile predicted by the 
linear response calculation is correct also outside the lin- 
ear regime (large Vo), but not all the way up to Vo — ► 00 
as on a sphere^, (b) This density profile n is also correct 
when Vo drops into the An(ro) oc V regime. However, it 
may be masked by density modulation due to finite size 
effects in small systems^. 

On a torus at at v = 1/3, the CM wavefunctions 
span a three-dimensional spaced. States in a homoge- 
neous system can be factorized into a CM and relative 
parts, * = ^cMiprei- Even the incompressible Laugh- 
lin ground state is thus triply degenerate on the torus 
and the electron density in this state strongly depends 
on which linear combination we choose for its ^cm- Im- 
purities lift the degeneracy but the splitting AE remains 
very small, AE <C Vo < E g . Density responses in each 
state Anv a — ny — ny- _>o arc n °t identical but differ- 
ences and vanish with increasing system size. The roles 
of the three states can be interchanged by moving r , 
Eq. [1] We always chose the state with lowest energy for 
plots in this article. 

We attribute the An(r ) oc V 2 behaviour to the situa- 
tion where the CM part of the wavefunction enters as a 
degree of freedom. Matrix elements of the electron den- 
sity between two states are proportional to the overlap of 
their CM parts. In particular, if these are mutually or- 
thogonal, the corrections to the density linear in Vo will 
vanish. For stronger impurity potentials, the CM degree 
of freedom may be frozen out and linear admixtures to 
the ground state wavefunction imply linear corrections 
to the electron density. For electrons on a sphere or im- 
purities discussed in the next section, the CM degree of 
freedom is never enabled and the density response is lin- 
ear down to Vq — > 0. 



II. RESPONSE TO A 5-LINE IMPURITY 

Comparing results for different system sizes is essen- 
tial in finite system studies. To keep such calculations 
tractable we must return to impurities which preserve 
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FIG. 3: The Laughlin state responding to a 5-line repulsive 
impurity. Responses in systems of different sizes are compared 
among each other and also to the linear response modelii for 
point impurities, cf. Fig. \T\ 
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FIG. 4: Response of the polarized | and | (almost identical 
to that of 3) states to a 5-line EI. Eight electrons considered, 
Vo = 0.01e 2 /(47T£^ ). 



some symmetries of the system. The 5-line impurity 

V s _\i nc {x, y) = V Q S(x - x ) , 

in a rectangle with periodic boundary conditions has a 
similar position as a point impurity on a sphere. The for- 
mer conserves translational symmetry along y, the latter 
conserves projection of the angular momentum to one 
axis. A straightforward consequence is that a 5-line im- 
purity does not mix the center-of-mass degenerate ground 
states on a torus. From the viewpoint of perturbation 
theory, the ground states are then non-degenerate and 
the regime of An(ro) oc Vq 2 , Fig. f2J is thus absent. Con- 
sequently, An(ro) is linear in Vq even for \Vo\ < E g . 

Throughout the rest of the article we will only use the 
5-line impurity with a strength of Vo = 0.01e 2 /(47re£o) 
which is for all studied states <C E g . Calculated re- 
sponses then provide information on the relative robust- 
ness of the different states (note however that the re- 
sponse is not a property of the ground state alone but 
depends also on low excited states). The form of the 
response is found to be independent on Vo as long as 
| Vo| < E g . Potentially, V^_j^ ne may be useful in studies 
of domain walls but we do not follow this aim in this 
article^. 

Turning to the -| Laughlin state, a linear-response 
analysis as in RefJ^ yields An(x,y) oc exp(— x 2 /s 2 ) for 
xq = 0. It is noteworthy that the form of An(x, y) should 
be now governed by the width s of the peak of x(q) rather 
than by its position qo. However, this simple model of 
x(q) seems to fail as the Laughlin state density response 
is again oscillatory, Fig. [3l yet with different period than 
for point impurities (dotted line). 

Keeping in mind the spin degree of freedom, an impu- 
rity can principally be one of the following three different 
types (5 CT | is a projector to states with spin |) 

N * { <W + <W ' Er 

H imp = E imp y2V(fi)(»i S^-8 ail 'MP . (2) 
<=i [ <5^ T 'DP' 



The electric potential impurity (EI) stands for a simple 
potential modulation (<5 CTj -f +^o-ii = 1)j the magnetic im- 
purity (MI) can be viewed as a spatially varying Zeeman 
field [V(fi) (X> (Scr^ — Scr z i) oc B z (fi)af] and the delta- 
plus impurity (DP) acts as a potential modulation seen 
only by spin-up electrons. The last type of impurity, DP, 
is less likely to occur in physical systems, however, it is 
helpful to understand the mechanisms governing inner 
structure of the states in study. 

It is important to check whether the response does not 
decay with increasing system size. Such a behaviour may 
imply that the response vanishes in large enough systems. 
However, in most of cases we find even a slight increase 
of the density response An(r), e.g. Fig. Exceptions 
from this will be explicitly stated. 



A. Electric potential impurity (EI) 

The polarized states of ^ , | and | respond all very sim- 
ilarly to a weak potential impurity. The singlet states at 
I and I respond more weakly and this is particularly 
apparent in the latter case. We give details on this ob- 
servation in this Section and then focus on the | singlet 
state in the succeeding Sections. 

It is not surprising that the Laughlin state and the 
polarized state at | have a similar response. They are 
particle hole conjugates (y and 1—v) and even though in- 
homogeneities break the symmetry of the Hamiltonian 14 , 
the effect appears extremely weak on the scale of Fig. fU 
Also the I polarized state and the Laughlin state show 
almost the same response, Fig. @] This applies both to 
its strength and the position of the first maximum (or 
node). Such a finding is non-trivial as the two states cor- 
respond to different filling factors in the CF picture. It is 
also non-trivial from the analogy between | and | . Here 
both states correspond to v* = 2 but the wavefunctions 
show different correlations on the electronic level^. 

Now turn to the singlet GS, Fig.0 The | singlet seems 



(a) I singlet (b) | singlet 

FIG. 5: The singlet states subject to a potential impurity. 
Response of the polarized states (thin dashed lines) is shown 
for comparison. 



to respond less strongly compared to the polarized state 
but the difference is small, Fig. [5p (crosses and the thin 
line). The difference between the polarized and singlet 
state at |, however, is striking, Fig. [3^- Measured by 
An(ro), the singlet appears about four times more robust 
than the polarized state. 

These findings can neither be completely explained 
with the intuition of non-interacting CFs nor purely on 
the basis of the numerically calculated gaps. For the for- 
mer, we would expect both the | and | singlet to have 
the same response as the Laughlin state. In all three 
states all CFs reside in the lowest CF Landau level and 
they are excited to the first CF LL by the impurity. On 
the other hand, the CF intuition is correct for the | and 
| polarized states. 

Regarding the gap energies it is also unexpected that 
the g spin singlet seems more robust than the polarized 
state whose gap is larger. An explanation for the stability 
of the singlet state must therefore lie in the structure of 
the wavefunctions rather than in the spectrum. 



III. SPIN PAIRING IN THE § SINGLET STATE 

The density-density correlation functions provide an 
alternative view at the internal structure of a many body 
state 1^). With spin degree of freedom, there are three 
distinct types denoted by g-ff (r), ga(r) and fffxM- With 
8 a *\ being a projector on spin up single-electron states, 
the first one is defined as (^\S(r — r' — r")6 <7 '^6 (T "^\'^}, 
the others analogously. 

The v = 3 singlet state displays remarkable struc- 
tures in these spin-resolved correlations, Fig. [51 While 
<7l"C (r) has a pronounced shoulder around r w 2£rj, the un- 
like spins have a strong correlation maximum near 3.5^o- 
Even though not shown in Fig. [6j this was checked for 
several different system sizes and the above dimensions 
turned out to be always the same. Note that these struc- 
tures are significantly stronger than in the cas e 14 ' 22 of the 
Laughlin state at v = ^ . We interpret this as a signature 
of spin pairing: two electrons with opposite spin form an 
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FIG. 6: Correlations between like spins and unlike spins in 
the | singlet state. Ten electrons in a square elementary cell 
(C marks its centre), g(r) along its diagonal is shown. Note 
the maximum in gn(r) around 3A£q. 



object of the characteristic size of approximatclly three 
magnetic lengths which is the mean intcrparticle separa- 
tion ^/2tt /vIq « 3.1fo- Qualitatively the same behaviour 
was found in the | singlet stated. 

Moreover, combining the correlation functions above 
into a spin-unresolved one 2 ^, g*(r) = g-\-\ + g-\i, we arrive 
at a result strongly reminiscent of the completely filled 
lowest Landau level (i.e. the ground state at v = 1), 
Fig. [7] Such a state is characteristic by a simple cor- 
relation hole at r = passing over monotonously to a 
constant^: g u= \(r) = 1 — cxp(— r 2 /2£%). Surprisingly, 
however, the correlation hole found in g*(r) corresponds 
to <jv = i with £\ replaced by 2£q with a rather good pre- 
cision (the fit in Fig. [7]). 

In a finite system with N e electrons, g* — g v —\ would 
mean that these electrons form the standard v = 1 state 
and behave as if they felt N m = N e magnetic flux quanta. 
With the replacement £\ — > 2£q we conclude that the elec- 
trons feel only N m = N e /2 flux quanta, or that each pair 
of electrons feels one flux quantum. Based on observa- 
tions in Fig. 16171 we suggest that the electrons in the | 
singlet state form pairs of total spin zero with charac- 
teristic size of 3^o and these pairs condense into a state 
resembling the completely filled lowest Landau level. 

The observed robustness of the singlet states against 
the potential impurities thus may be related to the ro- 
bust incomprcssibility of the filled lowest Landau level, 
possibly helped also by the fact that the relevant par- 
ticles are not single electrons but rather electron pairs. 
We now proceed with the investigation of the spin singlet 
states keeping in mind this picture providing a guidance 
at the interpretation. 



A. Magnetic impurity (MI) 

A magnetic impurity creates a relatively strong spin 
polarization p(x) = [n-f (a;) ~ n i( x )]/ n ( x ) of the singlet 
states, Fig. [BJ;,d. The density of electrons with spin f 
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FIG. 7: Data of Fig. [6] combined into a spin-unresolved cor- 
relation function of the | singlet state closely resemble the 
completely filled lowest Landau level with N e /2 particles (the 
fit). The g t (r) in the whole elementary cell and a section 
along the diagonal are shown. 



at x = decreases by 40% and 30% for | and |, re- 
spectively. On the other hand, a potential impurity of 
the same strength Vo changes the total density of these 
states only by 5% and 1.5%, Fig.O 

This observation is compatible with the concept of sin- 
glet pairs forming a v — 1 state (| singlet). While the 
potential impurity forces the whole unwieldy pairs to re- 
arrange, it is relatively easy to polarize the pairs without 
moving their centre of mass. In terms of the many-body 
states, the weak density response suggests that the cor- 
relations in the ground state are similar to those of low 
lying excited states. 

Several remarks should be added, (a) Magnetic impu- 
rity changes not only the polarization but also the den- 
sity, Fig. (5^,b. This effect is proportional to Vq and it is 
a consequence of that the density and spin density oper- 
ators do not commute within the lowest Landau level^s. 
Regardless of the sign of Vo, the density always increases 
at the impurity position, (b) The polarization profiles 
show only one node in our finite systems, Fig.[5V,,d. Char- 
acteristic length related to screening of magnetic impu- 
rities, if present at all, is therefore probably rather large, 
(c) The | singlet state with a magnetic impurity is the 
only case in this work where increasing the system size 
leads to a notably smaller density (and also polarization) 
response. 
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FIG. 8: Density (a,b) and polarization (c,d) response of the 
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2 and I singlet states to a magnetic impurity. 



to a potential and to a magnetic impurity even though 
EI+MI=2xDP on the level of the Hamiltonian, Eq. H 
We find that the density responses to a potential impu- 
rity (EI) and the DP impurity are in ratio ~ 4 : 1, Fig. [J^t 
while the polarization responses to the magnetic impu- 
rity (MI) and the DP are ps 3 : 1, Fig. [9]b. The both 
being more than 2 : 1 means that the capability of the 
individual spin species to answer to external perturba- 
tions is suppressed, hence suggesting that their mutual 
correlation is strong. It is also worth of recalling that the 
response of the | singlet state to a DP is thus more than 
an order of magnitude smaller than the one of the Laugh- 
lin state, Fig. A naive interpretation of v* = 2, for 
example, as of two independent v* = 1 systems (corre- 
sponding to Laughlin wavefunctions) with different spins 
therefore fails utterly. 



IV. SUMMARY 



B. Delta plus impurity (DP) 

Finally we investigate the | singlet state subject to the 
'delta plus' impurity, Eq. [2J and found weaker responses 
than for other types of impurities, Fig. [9] Put into the 
context of our singlet-pair incompressible state this find- 
ing reflects the strong correlations within the pair. 

We point out that the response to a DP impurity can- 
not be derived only from the knowledge of responses 



We find that the polarized v 



3 , | and | states all re- 
spond similarly to isolated impurities, represented in this 
paper by a weak 5-line impurity. Since | gives slightly 
different results than the other two systems we conclude 
that the particle-hole conjugation (| and |) is a stronger 
link than the reversal of the effective field in the compos- 
ite fermion picture (| and |). 

The singlet states react differently: both | and | 
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FIG. 9: Density (a) and polarization (b) response to all three 
types (Eq. [2J) of 5-line impurities acting on the | singlet 
(N e = 8). Density response of the Laughlin state to a po- 
tential impurity is shown for comparison in (a) by the thin 
line. 



respond quite unequally and more weakly to an elec- 
tric potential impurity than -|. In particular, | gives a 
much weaker response than |. This was unexpected be- 
cause the Laughlin state has the largest incompressibil- 
ity gap. The spin-resolved and spin-unresolved density- 
density correlation functions of the | singlet state suggest 
that electrons in it appear in zero-spin pairs with char- 
acteristic size of 3 magnetic lengths and these form a 
full-Landau-level-like state. High polarizability by mag- 
netic impurities and a relatively small effect of impurities 
affecting only one spin species are compatible with this 
interpretation. 
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